Quantum computing via defect states in two-dimensional anti-dot lattices 
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We propose a new structure suitable for quantum computing in a solid state environment: designed 
defect states in antidot lattices superimposed on a two-dimensional electron gas at a semiconductor 
heterostructure. State manipulation can be obtained with gate control. Model calculations indicate 
that it is feasible to fabricate structures whose energy level structure is robust against thermal 
dephasing. 



At present an intensive search is taking place for solid- 
state structures which are suitable for quantum comput- 
ing; a typical example consists of gate-defined double-dot 
systems studied by several groups [lj, |3, |3J, U, la, |6| . A 
necessary requirement for a practical application is scal- 
ability Q, and many of the existing structures do not 
immediately offer this possibility. Here we propose an 
alternative scheme: quantum-mechanical bound states 
which form at defects in an anti-dot superlattice defined 
on a semiconductor heterostructure. Scalability is not a 
critical issue for the suggested structures, which enable 
the fabrication of a large number of solid-state qubits 
with no particular extra effort. The flexibility offered by 
e-bcam or local oxidation techniques allows the sample 
designer to optimize the samples for many different pur- 
poses with a very high degree of control. 

Anti-dot lattices on semiconductor heterostructures 
have been a topic of intense research due to their inter- 
esting transport properties. In the semiclassical regime 
novel oscillatory features in magnetoresistance have been 
discovered @, and as the lattice spacing is diminished 
and the quantum regime is approached, exotic energy 
spectra, such as the Hofstadter butterfly [£| may become 
experimentally accessible. The fabrication of anti-dot 
lattices with lattice constants as small as 75 nm has been 
demonstrated in experiments |Toj . Smaller lattice con- 
stants are however expected to be within experimental 
reach [Tl| leading to a further enhancement of quantum 
effects. We shall in this paper demonstrate that state- 
of-the-art anti-dot lattices may have important practical 
applications in quantum information processing. 

Consider a two-dimensional electron gas (2DEG) at 
a GaAs heterostructure |l2j superimposed with a trian- 
gular lattice of anti-dots with lattice constant A. In the 
effective-mass approximation the two-dimensional single- 
electron Schrodinger equation reads 



2m 



r V* + £>(r-R< 



^n(r) = E n ip n (r), (1) 



where the sum runs over all anti-dots i, positioned at 



Each anti-dot is modelled as a circular potential barrier 
of height Vb and diameter d, i.e. V(r) = Vq for r < d/2, 
and zero elsewhere. It is convenient to express all energies 
in terms of the length scale A. Assuming that Vo is so 
large that the eigcnfunctions ip n do not penetrate into the 
anti-dots, i.e. ip n — in the anti-dots, Eq. Q simplifies 
to 



-A 2 V^„(r) =e n ip n (r), 



(2) 
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where we have introduced the dimensionless eigenener- 
gies e n = E n A 2 2m*/h 2 . For GaAs h 2 /2m* ~ 0.6 cVnm 2 . 

We first consider the perfectly periodic structure de- 
fined by the Wigner-Seitz cell shown in the left inset of 
Fig. n For definitencss, we now take d/A = 0.5. Im- 
posing periodic boundary conditions leaves us with the 
problem of solving Eq. @ on a finite-size domain. This 
class of problems is well-suited for finite-element calcula- 
tions, and the available software packages make the re- 
quired computations simple, convenient, and fast |14| . 
Fig. shows finite-element calculations of the bandstruc- 
ture along the high-symmetry axes indicated in the right 
inset of the figure. For state-of-the-art samples A ~ 75 
nm, implying a band-splitting of the order of 3 meV be- 
tween the two lowest bands at the T-point. On the figure 
we have also indicated the gap i? e g below which no states 
exist for the periodic structure. 

Next, we turn to the case where a single anti-dot has 
been left out of the lattice. Relying on the analogy with 
photonic crystal fibres, where similar ideas have been 
used to design confined electromagnetic waves [l^, we 
expect one or several localized states to form at the loca- 
tion of the 'defect'. The eigcnfunctions ip n corresponding 
to localized states decay to zero within a finite distance 
from the defect, and it is again sufficient to solve Eq. 

on a finite-size domain. The inset in Fig. [21 shows 
finite-element calculations of eigenfunctions correspond- 
ing to the two lowest eigenvalues for the geometrical ratio 
d/A = 0.5. The computed energy eigenvalues are con- 
verged with respect to an increase of the size of the do- 
main on which Eq. (0 is solved. The two lowest eigenval- 
ues correspond to localized states, whereas higher eigen- 
values correspond to delocalized states (not shown). The 
second lowest eigenvalue is two-fold degenerate, and we 
only show one of the corresponding eigenstates. One ob- 
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FIG. 1: (Color online) Bandstructure for the periodic struc- 
ture. The ratio between the diameter of the anti-dots and 
the lattice constant is d/A = 0.5. Only the five lowest bands 
are shown. On the (dimensionless) energy axis we have in- 
dicated the gap which can be considered as the height 
of an effective potential (see text). Left inset: Wigner-Seitz 
cell (grey area) for the periodic structure. Circles indicate 
anti-dots. Right inset: First Brillouin zone (grey area) with 
indications of the three high-symmetry axes along which the 
bandstructure was calculated. 



serves that the shown eigenstate does not exhibit the un- 
derlying six-fold rotational symmetry of the lattice. This 
can be traced back to the fact that the mesh on which 
Eq. (J2J was solved also lacked this symmetry. However, 
as recently shown by Mortensen et al. even weak dis- 
order in the lattice leads to a significant deformation of 
the higher-order eigenstates, and the shown eigenstate is 
thus likely to bear a closer resemblance to the states oc- 
curring in experimental structures, rather than the one 
found for an ideal lattice. Similarly, we note that the 
formation of defect states does not rely crucially on per- 
fect periodicity of the anti-dot lattice, which thus allows 
for a certain tolerance in the fabrication of the anti-dot 
lattice. 

Fig. |21 also shows finite-element calculations of the 
lowest eigenvalues corresponding to localized states as 
a function of the geometrical ratio d/A. In addition, 
the gap $ e ff as indicated on Fig. Q is plotted as a func- 
tion of d/A. The gap gives an upper limit to the ex- 
istence of bounds states and can be considered as the 
height of an effective two-dimensional spherical potential 
well in which the localized states reside. For GaAs with 
d/A = 0.5 and A = 75 nm the energy splitting of the 
two levels is AE = E^ — E\ ~ 1.1 meV which is much 
larger than fc^T at sub-Kelvin temperatures. Thus, a 
missing single anti-dot in the lattice leads to the forma- 
tion of a quantum dot with two levels at the location 
of the defect with an energy level structure suitable for 
a charge (orbital) qubit. As d/A is increased the con- 
finement becomes stronger and the eigenvalues and their 
relative separations increase. Moreover, the number of 
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FIG. 2: (Color online) Energy spectrum for a single quantum 
dot. The three lowest dimensionless eigenvalues, ei, £2, £3, 
(corresponding to localized states) as a function of the ratio 
between the anti-dot diameter d and the lattice constant A. 
The full line indicates the height # e g of the effective potential 
giving an upper limit to the existence of bound states (see 
text). The thin dotted line is the semi-analytic expression 
given in Eq. 0. Inset: Localized eigenfunctions ^1 (r) (upper 
panel) and ip2(jc) corresponding to the eigenvalues £1 and £2, 
respectively, for d/A = 0.5. The absolute square |-0* ( r ) 1 2 5 * = 
1, 2, is shown. 



levels in the quantum dot can be controlled by adjust- 
ing d/A, allowing for n = 1, 2, 3, . . . levels in the quan- 
tum dot. In particular, for any d/A < 0.42 a single-level 
quantum dot is formed. 



For sample optimizing purposes it is convenient to have 
simple expressions for the eigenvalues. In the limit of d/A 
approaching 1, the problem can be approximated with 
that of a two-dimensional spherical infinite potential well 
with radius A — d/2. For this problem the lowest eigen- 
value is e^ oo) = A 2 ag il /(A - d/2) 2 , where a ,i ^ 2.405 
is the first zero of the zeroth order Bessel function. Al- 
though this expression yields the correct scaling with d, 
the approximation obviously breaks down for small val- 
ues of d/A. In that limit we follow the ideas of Glazman et 
al. p"?| who studied quantum conductance through nar- 
row constrictions. The effective one-dimensional energy 
barrier for transmission through two neighboring anti- 
dots has a maximum value of it 2 , and we thus approxi- 
mate the problem with that of a two-dimensional spher- 
ical potential well of height tt 2 and radius A. The lowest 

(7T 2 ) 

eigenvalue e\ for this problem can be determined nu- 
merically, and we find e\ ~ 3.221. Correcting for the 



1ow-g?/A behavior we find 
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In Fig.[2]we show this expression together with the results 
for the lowest eigenvalue determined by finite element 
calculations. As can be seen on the figure, the expression 
given above captures to a very high degree the results 
obtained from finite-element calculations. For the higher- 
order eigenvalues similar expressions can be found. 

The leakage (transmission probability for penetrating 
the effective potential) due to a finite size of the anti- 
dot lattice can be found in the WKB approximation . 
Multiplying by a characteristic attempt frequency we get 
the following estimate for the inverse life-time 
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Td(E) 



E -2NA^ 



-(V oa -E) 



2m* A 2 



(4) 



where N is the number of rings of anti-dots surrounding 
the defect, and V e s = ^) e «h 2 /2m* A 2 . For GaAs with 
A = 75 nm, d/A = 0.4, and N = 1,2, 3, 4, 5, respectively, 
we find Td ~ 0.8 ns, 0.3 fis, 90 [is, 30 ms, 10 s. We 
see that even relatively small 'superlattices' offer nearly 
perfect confinement. 

We next consider the case where an anti-dot and one 
of its next-nearest neighbors have been left out of the 
lattice. Due to the close proximity of the resulting quan- 
tum dots, the different states of the two quantum dots 
couple with a coupling determined by the overlap of 
the corresponding single-dot wavefunctions. In partic- 
ular, for two single-level quantum dots, L and R, with 
corresponding states \L) and \R), respectively, a bond- 
ing |— ) = (\L) — \R))/y/2 and an anti-bonding state 
|+) = (|L) + \R))/\J~2 form. The corresponding eigenen- 
ergies are E± = E ± \t\ with E being the eigenenergy 
corresponding to each of the states \L) and \R), and t 
being the tunnel matrix element. From the eigenenergy 
splitting we easily obtain the tunnel matrix element as 
\t\ = {E + -EJ)/2. 

The coupling of the two levels can be tuned using a 
metallic split gate defined on top of the 2DEG in order 
to control the opening connecting the two quantum dots. 
By increasing the applied gate voltage one squeezes the 
opening, thereby decreasing the overlap of the two states 
\L) and \R). In the following we model the split gate with 
an infinite potential barrier shaped as shown on the inset 
in Fig. |2| Changing the applied gate voltage effectively 
leads to a change of the width w of the opening, which 
we in the following take as a control parameter. 

In Fig. |31 we show finite-element calculations of the 
dimensionless tunnel matrix element |t| = \t\A 2 2m* /ti 2 
as a function of the geometrical ratio w/A for a number 
of different values of d/A in the single-level regime, i.e. 
d/A < 0.42. For GaAs with A = 75 nm and d/A = 0.4, 
w/A = 0.6, the tunnel matrix element is \t\ — 0.015 meV. 




FIG. 3: (Color online) Coupling between two single-level 
quantum dots. The dimensionless tunnel matrix element |r| 
as a function of the ratio between the width w of the opening 
defined by the split gates and the lattice constant A for dif- 
ferent values of d/A (0.2, 0.3, 0.4) in the single-level regime. 
The width w is defined as the shortest distance between the 
split gates. Inset: Time propagation of an electron initially 
prepared in the state \L) (uppermost panel). Parameters are 
d/A = 0.4 and w/A = 0.6 which for GaAs with A = 75 nm 
implies as oscillation period of T = 0.14 ns (see text). The fol- 
lowing panels show the state of the electron after a time span 
of T/8,2T/8, 3T/8 (lowest panel), respectively. The absolute 
square |i/>(r)| 2 of the electron wavefunction is shown. 



With this coupling an electron initially prepared in the 
state \L) is expected to oscillate coherently between \L) 
and \R) with a period of T = h/2\t\ = 0.14 ns. We note 
that the period agrees well with the time scale set by the 
life-time obtained from Eq. (J2J with N = 1. According 
to the figure the coupling varies over several orders of 
magnitude, thus clearly indicating that the coupling of 
the two quantum dots can be controlled via the applied 
gate voltage. 

We have performed a numerical time propagation of 
an electron initially prepared in the state \L). In the 
inset of Fig. [3] we show a number of snapshots at differ- 
ent points in time as the electron propagates from the 
left to the right quantum dot. Once located in the right 
quantum dot, the electron starts propagating back to the 
left quantum dot (not shown), confirming the expected 
oscillatory behavior. 

Considering the double-dot as a charge qubit, one- 
qubit operations may be performed by controlling the 
tunnel matrix element as described above. Alternatively, 
one may consider the spin of two electrons, each localized 
on one of the quantum dots, as qubits. In that case the 
qubits (the spins) couple due to the exchange coupling, 
which again depends on the amplitude for tunneling be- 
tween the two quantum dots. In this manner one may 
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perform two-qubit operations as originally proposed in 
Ref. Q. 

In this work we have carried out a number of model cal- 
culations showing that an implementation of qubits using 
defect states in an anti-dot lattice is feasible. While we 
have here only considered the most basic building blocks 
of a quantum computer, a single charge qubit or two 
spin-qubits, we believe that the suggested structure can 
readily be scaled to a larger number of qubits. It is not 
difficult to imagine large architectures consisting of an 
anti-dot lattice with several coupled defect states and/or 
linear arrays of defect states constituting quantum chan- 
nels along which coherent and controllable transport of 
electrons can take place ^jj. We believe that the sug- 



gested structure, when compared to conventional gate- 
defined quantum dots, has the advantage that less wiring 
is needed. The individual antidots need not be electri- 
cally contacted, which in the case of conventional gate- 
defined structures may be a critical issue for large struc- 
tures consisting of many quantum dots. 

In conclusion, we have suggested a new structure which 
seems to offer many attractive features in terms of flexi- 
bility, scalability, and operation in the pursuit of achiev- 
ing solid state quantum computation. 
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ration of the manuscript, and T. Markussen for sharing 
his numerical codes with them. 
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